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The spin structure functions of the system of quasifree fermions on mass shell are studied in 
a consistently covariant approach. Comparison with the basic formulas following from the quark- 
parton model reveals the importance of the fermion motion inside the target for the correct evaluation 
of the spin structure functions. In particular it is shown, that regarding the moment Y\, both the 
approaches are equivalent for the static fermions, but differ by the factor 1/3 in the limit of massles 
, fermions (m <C po, in target rest frame). Some other sum rules are discussed as well. 

O . 

CM ■ 

C ' I- INTRODUCTION 

Ctf ■ 
Hi ' 

. ■ Measuring of the nucleon spin structure functions represents an important tool not only for better understanding 
of the nucleon internal structure in the language of the QCD, but also for better understanding of QCD itself. These 
. functions contain an information, which is a crucial complement to the structure functions obtained in the unpolarized 
J> ' deep inelastic scattering (DIS) experiments. 
lO . The polarized experiments are more complex and difficult than the unpolarized ones, nevertheless the last decade 
has brought remarkable results also for the nucleon spin functions from the experiments at CERN (EMC, SMC) and 
SLAC (E142, E143, E154, E155). And the new experiments are running (HERMES) or are being under preparation 
(COMPASS). The data on polarized pp collisions are expected from the collider RHIC. For the present status of the 
research in structure functions see e.g. [EJ, the overview and citation therein. The more formal aspects of the 
polarized DIS are explained in j| . 

Also the interpretation and understanding of polarized structure functions seem be more difficult. For an example, 
until now it is not well understood, why the integral of the proton spin structure function g\ is substantially less, than 
expected from very natural assumption, that the nucleon spin is generated by the valence quarks. Presently, there is 
Qh! a tendency to explain the missing part of the nucleon spin as a contribution of the gluons. It has been also suggested, 
that the quark orbital momentum can play some role as well Q- Q. 

The spin in general is a very delicate quantity, which requires correspondingly precise treatment. It has been argued, 
that for correct evaluation the quark contribution to the nucleon spin it is necessary to take properly into account 
the intrinsic quark motion - [Of. Necessity of the covariant formulation of the quark - parton model (QPM) for 
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the spin functions has been pointed out in 14 1. These requirements are not satisfied in the standard formulation of 
C$ ' the QPM, which is currently used for analysis and interpretation of the experimental data. 

In this paper we shall attempt to demonstrate the role of the intrinsic motion for the spin structure functions, using 
very simple model of the system quasifree fermions on mass shell. The basic requirement is consistently covariant 
formulation of the task for the system of fermions, which are not static, being characterized by some momenta 
distribution in the frame of their centre of mass. The spin structure functions of such system are obtained in Sec. || 



and the sum rules following from these functions are shown in Sec. HI. In the Sec. IV a comparison with the formulas 



of the standard QPM is done. The last section is devoted to the short summary. 

II. SPIN STRUCTURE FUNCTIONS IN COVARIANT APPROACH 

Let us imagine a system of three quasifree charged fermions with the spin 1/2 and mass m, for which the following 
conditions are satisfied: 
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1) The distribution of fermion momenta in the frame of their centre of mass is described by some spherically 
symmetric function G, 

J G( Po )d 3 p = 3; po = v/™ 2 +P 2 - (1) 
The free fermion states are described by the spinors 

^pA x ) = -7^ U (P' A)exp(-ipx); / tjA x (x)ip P! x(x)d 3 x = 1, (2) 



where f2 is the normalization volume and 

«(p,A) = -jL( J x . ); N=^^, 4^ = 1- (3) 

We assume 

^na<j)x = \<j)x, A = ±~, (4) 

which means, that the spin projection of the fermion in its rest frame is ±1/2 in the given direction n; |n| = 1. 
2) By G±x/2 we denote function, which measures probability, that fermion is in the state ^.±1/2: so that 

G( Po ) = G +1/2 (p ) + G_ 1/2 (po) (5) 

and we assume 

r AG(p )d 3 p=l; AG( Po ) = G +1/2 (po)-G^ 1/2 ip ). (6) 
The difference AG consists of the corresponding contributions Ahj from the three fermions: 

3 

AG(po) = ^2 Ah k(Po); Ah k (p ) = h ki+ i/ 2 (po) - h k _ 1/2 (po)- (7) 

k=l 

Later on, we shall need also the distribution 

3 

H(p Q ) = Y,etAh k ( Po ), (8) 
k=l 

where e k are the fermion charges. 

What is the resulting spin (total angular momentum) related to the whole system? Let us calculate the integral of 
the matrix elements 

(nj) = J jf £ Gx( Po ) (4>lx(x)*ti$pA*j) <? xd *P, (9) 
where the angular momentum j consists of the spin and orbital part 

]k = Sfc + h = ;r n _ - ISkhnPl-^ — • (10) 
2 V U <T k J dp m 

Since the total angular momentum j is a conserving quantity, which commutes with the term pa, a simple calculation 
gives 

ipl x (x)njip p ^x(x) = i (A + e k i m n k pix m ) . (11) 
So, after inserting to Eq. (j^) and using the assumption (^|) one gets 
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since the term £ki m nkPix m , due to spheric symmetry, vanishes. One can check, if n',n" are vectors, which together 
with n generate an orthonormal base in the frame of centre of mass of the three fermions, then a similar calculation 
gives 



Since we work with the probabilistic description (in terms of quantum mechanics with the statistical mixture of 
states) by means of the distributions G\, as a result we can obtain only the mean values of the total spin projections 
(J) = (0,0,1/2). Nevertheless one could consider a more rigorous (but more complicated) approach, in which the 
three fermion system is not constructed as the statistical mixture of plane waves, but as the composition of the 
three pure states j = 1/2, j z = ±1/2 with the condition, that the whole system represents a pure state J = 1/2, 
J z = 1/2. These states are represented by the relativistic spheric waves (spinors), which imply the corresponding 
probabilistic distributions G,G\,AG and H have spheric symmetry. In other words, if in our approach we assume 
the system in a pure state J = 1/2, then its probabilistic description in terms of the plane waves will be defined by 
the distributions G\, which are spherically symmetric. In fact, that is the reason, why we require spheric symmetry, 
deformed distributions G\ would contradict the eigenstate J = 1/2. 

Let us point out, in the relativistic case, having one fermion state with definite projection nj of the total angular 
momentum, one cannot separate its orbital and spin part (with exception of the special case when n || ±p), i.e. 
account with the fermion orbital momentum is crucial for a consistent calculation of the resulting spin. On the other 
hand, the similar calculation, in which the orbital part 1 is ignored, gives 



(n'j) = (n"j) = 0. 



(13) 



Obviously, the simplest way is to use the base like: 



n= (0,0,1) 



n' = (0,1,0) 



(1,0,0). 



(14) 






Since 
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(15) 



i=i 



one can write 
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where inserting the formula (p"5|), we take into account, that due to spheric symmetry the terms piPj (j ^ i) vanish 
and the terms p 2 can be substituted by p 2 /3. The last relation can be further simplified: 

<nS)= i|i GW (i + ^) d 3 P a. (16) 

One can observe, that the correspondence with Eq. (|l2|) takes place only for the system of static fermions. 

For further consideration, it will be useful to substitute the vector n, representing the direction of the fermion 
polarization, by the corresponding covariant polarization vector w CT (A), which satisfies 

w 2 {\) = -l 1 w(\)-p = (17) 

and 

«;(A) = A(0,n); \ = ±\ (18) 

in the fermion rest frame. The explicit representation of the vectors w(X) will be defined hereinafter. 

Now, let us expose this system as a (fixed) target to the beam of polarized electrons (e.g. helicity = +1/2) coming 
with the momentum 

k = ^ko,^k 2 -mlO,0^j (19) 

and let us calculate the form of corresponding differential cross-section. The spin dependent part of the cross-section 
for interaction with a single fermion in one photon approximation has the form 

d*~-L a f>W(q,8)T<f. (20) 
The antisymmetric tensor L a/3 ( A \ (see e.g. El) related to the electron beam reads: 

L a PM = m e s a p Xa s x q°, (21) 
where m e is the electron mass, s denotes its polarization vector 
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k 2 -m 2 e ,k , 0,0 ; s 2 = -l, ks = (22) 



and q = k — k' is the photon momentum. The antisymmetric tensor T a @( A ) related to the single fermion inside the 
target has a similar form: 

T a PW =me al3 x„q X w a {\), (23) 

where m and w(X) denote the fermion mass and polarization vector. If one assumes, that the electron scattering can 
be described as the incoherent sum of the interactions with the single plane waves, then the tensor T a ^^ reads 

T { a f = e apXa q x m [ H( Po )w°6((p + q) 2 - m 2 )^. (24) 
J Po 

Here the charge factors and the two possible signs of w 17 are included into the tensor through the distribution (||) . By 
the symbol w a we mean w a (X = +1/2). Let us remark, this form of the antisymmetric part of the hadronic tensor is 
very similar to that used in [ fl4"|| . Further, we can modify the 5— function term: 

6((p + qf ~ m 2 )d 3 p = 5(2pq + q 2 )d 3 p = + ^)d 3 p, (25) 

where £ is arbitrary constant, which only rescales the integration variable. Now, let us imagine, that our target is a 
part of the greater system, which is at rest with respect to the given reference frame and has the mass M, but at the 
same time the probing electron interact only with the three fermions. If we put 

£ = Mq = MV, (26) 



4 



then in the 5— function one can identify the terms known from the formalism of deep inelastic scattering: 



= (27) 

which is the Bjorken scaling variable, its value can be directly determined using only initial and final momenta of 
the scattered electron. This variable is in the S— function compensated by the ratio pq/Mv, which after boosting 
the whole target of mass M to the infinite momentum frame approximately represents ratio of dominating momenta 
components p' /P' of the fermion and the target. 

The explicit form of the polarization vector w can be found as follows. First, let us transform the vector w = (0, n) 
from the fermion rest frame to the target rest frame. After decomposition of the vector n to longitudinal and 
transversal parts with respect to the momentum fermion p, the corresponding Lorentz boost gives 

(0,n)^=fc n+ P^ p ). (28) 
\ m m(m + po) J 

Secondly, let us make a Lorentz boost of the whole target with mass M to some another frame, which is defined by 
the new components of the target momentum 

(M,0,Q,0)->P=(P 0) P); P 2 = M 2 . (29) 

Next, if we define the covariant vector S by its components in the target rest frame as 

5=(0,n), (30) 

then the polarization vector w can be written in manifestly covariant form 

w a = AP a + BS a + Cp a , (31) 

where A,B ) C are invariant functions (scalars) of the vectors P,S,p. These three functions are fixed by two the 
conditions (17) and by the constraint ( |28|) valid in the target rest frame. A simple calculation gives: 

A = - J S - , B = l, C=—A. (32) 
pP + mM m 



So, we have obtained explicit covariant form of the polarization vector w entering the tensor (|2J|), which can be now 
in accordance with the relations (p5|)-(p7|) rewritten 

rp(A) \ m f ( pP\ ax fpq \ d 3 p 



where we use the invariant term Pq instead of Mv and H(pP/M) instead of H(po) 



On the other hand, in accordance with the general rule (see e.g. 0]), the antisymmetric tensor appearing in 



the formula for the cross-section (EOT), has the form 



T ( a f = S a p Xa q X |m5 ct G! + [{Pq)S° - (qS)P^]^ , 



(34) 



where M, P, S represent the target mass, momentum and spin polarization vector, which satisfies 

S 2 = -1, PS = 0. (35) 

The invariants G\ and G2 are the spin structure functions. In the next we shall identify the parameters M, P, S in Eq. 
([m|) with those in the model described above and simultaneously we shall attempt to determine the spin structure 
functions corresponding to our target. First of all, we modify the Eq. ([mJ) by the substitution 

G s = MGi + §G 2 , G P = ^G 2 , (36) 
M M 

which gives 
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T 



(A) 



a/3 



^0Xaq X {S a Gs-P a G P }. 



(37) 



Comparison with Eq. ( J33] ) gives the equation for the structure functions: 

W (S°G S - P-G r} = / H (£) •« - ,) J. 

Because of the antisymmetry of the tensor £ and after inserting from the relation (^l]) it follows that 



S a G s - P a G P = 



2Pq 



(38) 



(39) 



where D is some scalar function and the functions A,B,C are given by the relations (|32|). After contracting with 
P a , S a and go- one gets the equations for unknown functions Gs, Gp and D: 



M 2 G P = 



2Pq 



Gs 



2PqJ \M) K P ' \Pq 



■ X 



d ^ + D.qS, 



Po 



(40) 



(41) 



qS-Gs~Pq- G P 



(42) 



Pq J po 



and inserting Gp, Gs from the first two equations to the last one gives the condition for D: 

m 



where we denote 



u = q + (qS) S - &-R 



(43) 



M 2 

Finally, inserting D from this equation to Eqs. (|o|), (|4l]) gives with the use of relations ( |32| ) the structure functions 

'pP\ pS 



G P 



2Pq 



H 



M J pP + mM 



1 + ^M\ PP -- Pq 
mM \ qu 



Pq 



Pi) 



(44) 



— / H (— 

2PqJ \M 



pS M / _ pu s 
pP + mM to V qu 



Pq 



Po 



(45) 



The s pin structure functions in the standard notation g\= M ■ Pq ■ G\, g% = (Pq) /M • G-i can be now obtained from 
Eqs. Mj: 



.'/i = r<i{G s -^G P 



<]2 



(Pq)' 
q s 



■G P , 9i+92 = PqG s , 



(46) 



where the functions Gs, Gp are given by relations (f44|), ([45|). Corresponding integrals, as shown in the Appendix, can 
be sim p lified to the form (A14), (A15). Let us remark, resulting functions 51,32, after inserting from the relations 
(A14), (A15) into Eq. (46) do not depend on the variable qS despite the fact, that such terms are present in the 
starting integrals (Q), (45) in a non-trivial way. This is a consequence of spheric symmetry of the distribution H, 
which as we have suggested, follows from the requirement J = 1/2. 
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III. SUM RULES 



For next analysis of the obtained structure functions it is convenient to express the integrals (|44|) , fl45f ) in the target 
rest frame, where P = (M, 0,0, 0) and S = (0,n). Detailed calculation is done in the Appendix. Now, let us assume 
Q 2 > 4M¥, then 

M = v/l + 4M 2 x 2 /Q 2 -> 1 



and using the second relation ( [46| ) and Eq. (A7) one gets 
r2 = / 92{x)dx = -it H(p ) I pi + 



Pi-Pt/ 2 \ g ( Po +Pi _ \ PrdpxdpT 



Po + m J \ M 



(47) 



n / H (Po) Pi 



Pi -Pt/ 2 \ Prdpidpr 
Po+m J po 



In the last integral, due to spheric symmetry of the distribution H, the terms proportional to p\ and p\ —p^/2 vanish, 
insofar that 

r 2 = 0, (48) 
which is the known Burkhardt-Cottingham sum rule . Similarly the third relation ( f46| ) and Eq. ( |A8| ) give 

Pt \ Prdpidpr 



(gi(x) + 92{x)) dx = ir \ H(p ) m + 



2 (p + m)J po 



After the substitution 



2irpxdp\dpT = dp 



and using relation ( f48|) one gets 



Ti = / 9i{x)dx = - / H(po) m + — ■ — 

2 ./ V 3 (po + m) ) po 



d 3 p 



Simple modification then gives 



ri = \ J h(po) 



1 2m 
3 3p 



d 3 p. 



(49) 



(50) 



(51) 



(52) 



More detailed analysis of t his r esult will be done in the next section. 

The relations (A7) and (A8) can be used also for the calculation of the higher momenta. Generally, if F is a 
function defined as 



F(x) = J K(p)S (*±* - d*p, 



then 



J x n F{x)dx = J J K(p)x n S 



^° ^ ^ 1 — x ) d 3 pdx 
M 1 ' 



M 



Po+Pi\ § (Po+P 1 _ x]d s pdx 



M 



i^n^pA d*p. 
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Application of this rule to Eqs. (A7) and ( |A8| ) gives after the substitution ( |5C| ) and with the use of the second and 
third relation (ffHI): 

J xg 2 dx = --^ J H(p ) (p - ^-^J d 3 p, (53) 
J x {g x + gi) dx = / H( Po ) [p Q + 2m) d 3 p. (54) 



These equalities imply relation 

J x ( 9l + 2g 2 ) dx = -^ j H{ Pa ) (2m + ^ d 3 p, (55) 
which in the limit of the massless fermions coincides with the Efremov - Leader - Teryaev (ELT) sum rule |l6| ] : 

J x( gi + 2g 2 )dx = 0. (56) 



IV. DISCUSSION 

In the previous sections we have studied the properties of the spin structure functions related to the system of 
quasifree fermions on mass shell. This system can be compared with the naive QPM, which is with embedded QCD 
corrections yet the basic tool for the analysis and interpretation of polarized and unpolarized deep inelastic scattering 
data. What is the difference between our approach and the naive QPM, if one speaks about the proton spin structure 
functions? To simplify this discussion, let us assume: 

1) Spin contribution from the sea of quark-antiquark pairs and gluons can be neglected. Then the three fermions 
in our approach correspond to the three proton valence quarks. So, in this simplified scenario, the proton spin is 
generated only by the valence quarks. 

2) In an accordance with the non-relativistic SU(6) approach the spin contribution of individual valence terms is 
given as 

s u = 4/3, s d = -1/3. (57) 

Let us point out, in the given context the term valence quarks means nothing else, than the three fermions with 
defined momenta distribution, charge, mass and polarization. 
Then according to the naive SU(6) version of the QPM we have 



2\ 2 2 , , /1\ 2 1 



.'/!(•'') = -> ( : olio) 3 U «ol(aO - (^3 J ^dvai(x)J , (58) 

corresponding to two the quarks with distribution u va i (x) and the one with distribution d va i (2) , which are normalized 
as 



u va i(x)dx = d va i(x)dx = l. (59) 



It follows, that 



T 1 =f gi(x)dx = A = 0.28. (60) 

This number overestimates more than twice the experimental value. Disagreement is generally interpreted as a 
contradiction with the assumption, that the proton spin is generated only by spins of the valence quarks. 

Now let us calculate the Ti in our approach. Let us denote momenta distributions of the valence quarks in the 
target rest frame by symbols h u and hd with the normalization 
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\ J K(j>o)d 3 p = J h d (p )d 3 p = 1. (61) 
These distributions are connected with the u va i(x) and d va i{x) defined above by the relation 

*,«,(») = J h g (po)6 ( P °^^ l|ql - x\ d 3 p. (62) 
The charge weighted distribution (||), in an SU(6) picture, reads 

H (Po) = E e ' A/l i(Po) = ((jf ~Mpo) Q) 8 ^d(Po)J • (63) 
Now, for simplicity let us assume the same shape of the distributions for both the flavours: 

^K(po) = h d (p ) = h(p ). (64) 



Then it follows 



G( Po ) = 3fc(po), AG(p ) = Mpo), ff(po) = ^(Po) (65) 



and the relations (|l6|) and (S3) can be rewritten: 



(n S ) = l/, W (I + ^)A, (66) 

r -l/^»(MS) dV (67) 

These relations imply: 

a) Because the distribution h has the defined normalization, the corresponding integrals reach their maximum in 
the limit, when the fermions are static (po = to). On the other hand in the limit of massless fermions (to <C po) these 
integrals represent only one third of their maximal value. In particular, the Ti satisfies: 

£ > r, > |. (68) 

b) Both the integrals are (up to the factor 5/9) equal. It follows, that in the case of non static fermions the Ti 
"measures" o nly the contribution from their spins, which is only part of the their angular momenta, see derivation of 



the relation (16). Fermions with momentum p ^ 0. which is not parallel to ±n, necessarily contribute to the total 
angular momentum also by some orbital part. 
Further let us notice, if we denote 

1elt = J x(gi + 2g 2 )dx, (69) 



then Eq. ( p5[) and the third relation QS5J) imply 

2 18 

m < — lELT ■ M < to. (70) 

3 5 

Why these two very simple approaches for description of the target consisting of the three fermions differ so strongly 
regarding the prediction Ti? The reason is following. The standard formulation of the QPM is closely connected with 
the preferred reference system - infinite momentum frame (IMF). The basic relations between the distribution and 
structure functions like 



9i(x) = ^^2e 2 j Aq j (x), F 2 {x) = x^e\qi{x) (71) 
are derived with the use of approximation 
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Pa 



= xP c 



(72) 



which seems to be plausible in the IMF. Nevertheless, in the covariant formulation this relation is equivalent to the 
assumption, that the quarks are static with respect to the proton, since the velocities Pj/po and Pj/Po are the same. 
In the proton rest frame it means p = 0. That is why both the approaches are equivalent for the static quarks but 
differ for the quarks, which have some intrinsic motion inside the proton. In our approach we do not use assumption 
((72|) and as a result if p a ^ xP a we obtain different relations between the distribution and structure functions. In 
other words, the fact, that the experimental value T\ is substantially under the value predicted by the naive QPM in 
standard formulation, can be in our approach interpreted as a direct consequence of the quark intrinsic motion. 

Of course, the approach discussed above concerns the simplified scenario of the quasifree fermions on mass shell. 
Naive QPM represents only a first approximation for a description of real nucleon, but the consistent accounting for 
the quark intrinsic motion as suggested in our approach can, in some aspects, improve this approximation considerably. 

Nevertheless, in the realistic case of partons inside the nucleon the situation is still much more delicate. The 
interaction among the quarks and gluons is very strong, partons themselves are mostly in some shortly living virtual 
states, is it possible to speak about their mass at all? Strictly speaking probably not. The mass in the exact sense 
is well defined only for free particles, whereas the partons are never free. However one can assume the following. 
The relations obtained in the previous sections can be used as a good approximation even for the interacting quarks, 
but provided that the term mass of quasifree partem is substituted by the term parton effective mass. By this term 
we mean the mass, which a free parton would have to have to interact with the probing photon equally as the real, 
bounded one. Intuitively, this mass should correlate to Q 2 : a lower Q 2 roughly means, that the photon "sees" the 
quark surrounded by some cloud of gluons and quark- antiquark pairs as a one particle - by which this photon is 
absorbed. And on contrary, the higher Q 2 should mediate interaction with more "isolated" quark. Moreover, one 
should accept that the value of the effective mass can even for a fixed Q 2 fluctuate. Such phenomenological model 
was suggested in [ [l3|| , but unfortunately calculation was based on the form of quark polarization vector which is 
not correct. Despite of that, the general considerations in mentioned paper can be sensible. Corresponding numeric 
recalculation with the correct input obtained in the present study for the invariants A,B,C,D [relations (j32|),([f3"l)1 
should be done in a separate paper. 



In the present paper we have studied the spin structure functions of the system of quasifree fermions on mass shell 
and with spherically symmetric distribution of their momenta. The main results can be summarized as follows: 

1) Using consistently covariant description of this simple system, we have shown how the structure functions depend 
on the intrinsic motion of the fermions. In particular, we have suggested, that the momenta T\ corresponding to the 
two extreme scenarios, of the static (massive) fermions and massless fermions, can differ significantly: Ti(m <C 



2) We have shown, what sum rules follow from the obtained spin structure functions. Further we have shown, how 
these rules are related to some sum rules well known from the QPM phenomenology. 

3) We have done a comparison with the corresponding relations for the structure functions following from the 
standard formulation of the naive QPM. Both the approaches are basically equivalent for the static quarks. Differences 
for quarks with intrinsic motion inside the proton are result of the conflict with the assumption p a = xP a , which is 
crucial for derivation of the relations between structure and distribution functions in the standard QPM. 

4) The difference between the experimental value T± for the proton and the corresponding value expected from the 
naive QPM, or at least a part of this difference, can be interpreted as a consequence of the quark motion inside the 
proton. 

Acknowledgement: I would like to thank Anatoli Efremov and Oleg Teryaev for many useful discussions and 
valuable comments. 



V. SUMMARY AND CONCLUSION 



Po)/Ii(po ~m) = 1/3. 



APPENDIX A: CALCULATION OF THE INTEGRALS RELATED TO G P ,G S 



Integrals in the relations (H), ( fig ) expressed in the target rest frame read 




(Al) 
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pn 



Po + m 



Po 



pq — (pn) |q| cosw 



q 2 sin 2 lj 



d 3 p 



\Mv J po 



2Mv 



H( PQ ) 



(A2) 



1 



Pn 1 ( pq - (pn) |q|cosw 

pn — -= |q| cosw 



q 2 sin 2 w 



Po + m to 

where cosoj = qn/ |q| . For integration we use the orthonormal system in which 



\Mv J Pq 



p = piei + p 2 e 2 +P3e 3 , ex 



q n-(nex)ex 

i — r, e? = — 1 e-? = ei x e?. 

|q|' ^-(nei) 2 



so one gets 



pq=-px|q|, pn = -pxcosw + p 2 smw, 



coscj 



qn 



(A3) 



(A4) 



After the substitution p 2 = Pt cos tp, p 3 — pt sin tp and taking into account that the terms proportional to cos ip 
disappear, the integrals can be rewritten 



cosw . v p\-p\ COS 2 if 

Gp = TTTTt / H(po) pi + 



2M 2 v 



|q| Po + m 



(A5) 



x g i Pqv+Pi |q| _ ^ p T dpidp T dtp 



Mv 



Po 



p 2 -, COS 2 </5 
TO (po + TO-) 



p a v + pi |q| 



prdpidpxdip 
Po 



(A6) 



where po = \/ m 2 + p 2 ^ + p 2 ■ Integration over (/? gives 



ttcosw f / zv p 2 -p|/2\ /po^+Pi q \ prdpidpr 

Gp = M2 / ^(Po Pi + -pr -7 — 5 T7 a; , 

M l v J V q Po+m / \ Mv J po 



(A7) 



G & 



Mv 



J H( Po ) (l 



p 2 T /2 



m (p + m) 



, po^+pi |q| _ ^ \ p T dpidp T 



Mv 



Po 



Further, using the relation 



M _ y J\ + AM 2 x 2 /Q 2 



(A8) 



(A9) 



one can check, that the argument of 8— function equals zero for 



Pi = Pi 



Mx — vrOp /Mx 



y/1 + 4m 2 T /Q 2 + v /l + 4Af 2 x 2 /Q 2 



2 _ 2,2 
m T = TO + p T . 



(A10) 



This is the first root of the corresponding quadratic equation, the second one is excluded, since in the effect of the 5- 
function this root is compatible only with negative energy p . The energy corresponding to the root (A10) is 



Po = Po = Mx 
Then in an accordance with the rule 



Pi q 



v 



= Mx - pi y/\ + AM 2 x 2 /Q 2 . 



(All) 
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5(f(x))dx = J2 %TzSr dx, f( Xj )=0 (A12) 



the 5— function in the integrals can be rewritten 



s /W + Pildl _ \ d = MSjpr-p^ 

1 Mu J p 1 /p Q + ^/l + AM 2 x 2 /Q 2 



and afterwards the integrals are simplified 

:cos ^ r Tmax F / 

Mi/ 7o ^ ■ | q | p Q + m J p 1+ p oy /i + 4i M 2 x 2 /Q 2 



G S = ™ f T ~ fffo) f i + - ) ^ (A15) 

^ io V m (po + m) y p x + p + 4M 2 x 2 /Q 2 

where pi and depend on pt according to Eqs. ( |A10| ) and ( |All ). For th e num eric calcu lation one should know the 
upper limit prmax for given x, Q 2 and po max . After inserting p\ from Eq. (A1C) into Eq. (All) one gets equation for 



PO max — M x Mx — jMx 



^l + AM 2 x 2 /Q 2 ^l + im^/Q 2 + ^/l + AM 2 x 2 /Q 2 



Instead of it is useful to solve this equation first for y = y/l + Arn^/Q 2 obtaining the two roots 



(A16) 



A ± y/A 2 + 4a(p max + a) Pomax-Ma: Q 2 

V+ = - 1 A = — ■ a = . A17 

2a ^/l + AM 2 x 2 /Q 2 4Mx 

Since y_ < 0, this root is excluded. The second root j/+ after some computation implies 

TO Tmax = Mx(2p 0max - Mx) + ^+Q2 j ^H 2 ' ?T m ** = \/ m T max ~ m2 ■ ( A18 ) 

In this way we have the recipe how to calculate the integrals related to the structure functions Gp, Gs corresponding 
to the distribution H(p )d 3 p. 
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